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Abstract 

We study the diffusion-driven kinetics of phase separation of a symmetric binary mixture 
(AB), confined in a thin-film geometry between two parallel walls. We consider cases where 
(a) both walls preferentially attract the same component (A), and (b) one wall attracts 
A and the other wall attracts B (with the same strength). We focus on the interplay of 
phase separation and wetting at the walls, which is referred to as surface- directed spinodal 
decomposition (SDSD). The formation of SDSD waves at the two surfaces, with wave- 
vectors oriented perpendicular to them, often results in a metastable layered state (also 
referred to as "stratified morphology"). This state is reminiscent of the situation where 
the thin film is still in the one-phase region but the surfaces are completely wet, and hence 
coated with thick wetting layers. This metastable state decays by spinodal fluctuations 
and crosses over to an asymptotic growth regime characterized by the lateral coarsening of 
pancake- like domains. These pancakes may or may not be coated by precursors of wetting 
layers. We use Langevin simulations to study this crossover and the growth kinetics in 
the asymptotic coarsening regime. 
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I. INTRODUCTION 



Consider a binary mixture (AB) with a miscibility gap, such that phase separation 
into A-rich and B-rich phases occurs below the critical temperature T c . If this 
mixture is quenched from a homogeneous state in the one-phase region into the two- 
phase region below the critical point, phase separation proceeds by the emergence 
and growth of regions enriched in either component. In the bulk, this process of 
spinodal decomposition or domain growth has been intensively studied |jl Q, El B, El ■ 
In recent years, the effect of surfaces on this behavior has become a focus of 



research 
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12, 



131 ] . It is often the case that a surface (S) has a 



preferential attraction for one of the components (say, A) of the mixture. In the 
one-phase region of the mixture, this attraction leads to the formation of surface 
enrichment layers of the preferred component. The thickness of these layers (which 
is comparable to the correlation length of concentration fluctuations in the bulk) 
becomes long-ranged at the critical point of the mixture jl^l . In addition, if the 
composition of the bulk mixture coincides with the B-rich branch of the coexistence 
curve describing phase separation in equilibrium, the surface is coated with an A- 
rich wettin g la yer for temperatures below the wetting transition temperature T w 

y y yy y 21, y . Of course, wetting ayers of mesoscopic thickness (in 

the idealized case, the thickness of the wetting layer diverges when complete wetting 



occurs 19, 2 
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) can only occur in macroscopic systems. Typic ally , theoretical 



approaches consider the idealized case of a semi-infinite geometry 



18, 



19 



20 



5 y, y, y, 



2l|. 



In a thin film of thickness D, finite-size effects significantly constrain the growth 
of wetting layers, e.g., for short-range forces between the walls and the A-atoms, the 
equilibrium thickness of a wetting layer is 0(\nD) y,y>E)l3- Further, in a thin- 
film geometry, the interplay between surface effects and bulk phase separation leads 



to a distortion of the phase diagram describing lateral phase separation 



the walls, analogous to the phenomenon of capillary condensation (22 , 
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28|,|29|. Moreover, the critical behavior changes its character from three-dimensional 
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to two-dimensional 30] . The interplay of wetting and phase separation in a thin film 
results in a rich phase behavior, and the equilibrium phase diagrams of thin fi 



have been investigated for a range of possible surface forces 
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. In the low-temperature region, relevant for deep quenches to a temperature 
T < T w , the typical state in a thin film with symmetric walls (i.e., both attract A 
with the same strength) is a laterally-segregated state with a B-rich pancake-like 
domain. This domain is circular with a macroscopic radius R, which is of the same 
order as the macroscopic linear dimension L parallel to the walls - see Fig. ^a). For 
a volume fraction 0b of B-atoms lying between the two branches of the coexistence 
curve 0sc OCX (-D), 0Bcoex(^) °f the thin film, the lever rule requires 



>B 



*)<f>%]o*JP) + 4 2 Lw , (i) 



where x is the volume fraction of the B-rich phase. For a film of volume L 2 D 
(Fig. ^) with periodic boundary conditions in the directions parallel to the walls, 
we have x = ttR 2 /L 2 . Assuming a strongly segregated state where 0Bc Oex (-D) — 0, 
^Bcocx(^) — 1, we obtain R ~ L^<f) B /ir for the radius of the B-rich domain. 

In this simple argument, we have assumed that the pancake is cylindrical and 
its surface is perpendicular to the walls. In fact, the interface is curved because (in 
the D — > 00 limit) the A-B interface meets the wall at a contact angle 9, given by 
Young's equation j^: 

a cos 6 = 7s - 7a • (2) 

Here, a is the interfacial tension between the A-rich and B-rich phases; and 7,4, 
7b are the surface tensions between the A-rich and B-rich phases and the wall, 
respectively. The surface is partially wet (PW) when a > 7b — 7a 
21, Q When T — > T w , 6 — > 0, and for T > T w one has a < 7b — 7^. In 
this complete wetting (CW) situation shown in Fig. HJb), there is no direct contact 
between the B-rich phase and the wall - the B-rich pancake is encapsulated by A- 
rich wetting layers for D — » 00. (For finite D, only precursors of wetting layers are 
possible.) The cross-sections of these states [Figs. [TJa)-(b)] are reminiscent of the 
states considered by Liu et al. |38| in the context of phase separation in cylindrical 



pores. Liu et al. made a distinction between elongated plugs in the PW case and 
capsules in the CW case, respectively. Note, however, that the lateral size of these 
domains in equilibrium for small D and L — > oo does not diverge with L, unlike 
the present d = 2 cases. Rather the lateral size t ~ exp[<T7r D 2 /(4k bT)], where D is 
the pore diameter. We also note that, in the CW case the thickness of the B-rich 
domain differs from D only by a small amount, that will be estimated below. 

It is also interesting to consider equilibrium morphologies in antisymmetric films, 
i.e., one wall attracts A and the other attracts B with equal strength. In Fig. ^c), 
we show a schematic of the PW state in an antisymmetric film. This state also 
resembles a pancake, except that the cross-section is trapezoidal in this case due to 
the difference in contact angles at the two walls. In Fig.^d), we show a schematic of 
the CW state in an antisymmetric film. In this state, there is a single A-B interface 
which is parallel to the walls. 

In this paper, we are interested in the kinetic processes which lead to the for- 
mation of these structures subsequent to a deep quench from the high-temperature 
disordered state. In these surface- directed spinodal decomposition (SDSD) processes, 
one typically observes first the formation of a layered structure, consisting of an A- 
rich layer followed by a depletion layer in A, etc. This layered profile propagates 
into the bulk, but at a later stage it may break up into a laterally inhomogeneous 
structure. The kinetics of such growth phenomena has important technological ap- 
plications, including the fabrication of nanoscale patterns and layered structures. 




However, despite much effort (for reviews, see |9|, [l0|, [ill, Q, 13() these phenomena 
are still incompletely understood. 

In the present paper, we will study the dynamics of these structure-formation 
processes for the case where hydrodynamic effects can be disregarded, so that the 
kinetics is purely diffusive. Therefore, our results will be applicable to solid mixtures, 
and the early stages of phase separation in polymer and fluid mixtures. (In a subse- 
quent paper, we will elucidate the role of hydrodynamic effects in these problems.) 
We shall consider the cases of both symmetric and antisymmetric films. This paper 
is organized as follows: In Sec. II, we describe our dynamical model for segregation 
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in confined geometries. In Sec. Ill, we focus on the kinetics of phase separation in a 
symmetric film, while Sec. IV studies segregation kinetics in an antisymmetric film. 
We conclude this paper with a summary and discussion in Sec. V. 



II. DYNAMICAL MODEL FOR SEGREGATION IN FILMS 



The model forphase separation at surfaces has been developed previously by 
Puri and Binder jl). In dimensionless units Q, 3], the free energy functional for 
an unstable binary mixture in a film extending from z = to z = D is j3| 



dr 



-T + T + I ( ^ )2 + V ^ 



+ 



dp{ -|[^(p,o)] 2 - M(£o) - T ^(p,o)^ 



Si 



+ 



J dpi -| W, D)f - D) + 7 ^(p, D)^ 



s 2 

Fb + Fs-l + F S2 



(3) 



Here, ip(r) is the order parameter which is proportional to the density difference 
between the two species, if)(r) oc Pa{^) — Pb{^)- It is normalized such that the coex- 
isting A-rich and B-rich bulk phases for T < T c correspond to ip — ±1, respectively. 
The first term on the RHS of Eq. (jHJ) is the bulk free energy Fb, with V(z) being the 
z-dependent potential due to the surfaces Si and 5*2. In our subsequent discussion, 
we will consider power-law potentials: V(z) = —Vq[(z + l)~ n ± (D + 1 — z)~ n ], 
where the + and — signs denote symmetric and antisymmetric films, respectively. 
Such potentials are common in the context of surface-molecule interactions, e.g., 
n = k — d with k = 6 and 7 corresponds to cases with non-retarded and retarded 
van der Waals' interactions, respectively j^j. The potentials are taken to originate 
behind the surfaces so as to avoid singularities at z = 0, D. 

The second term Fg x on the RHS of Eq. © is the local contribution from the 
surface Si located at z = 0. We have written f— (p, z), where p denotes the (d— 1) 
coordinates parallel to the surface, and z denotes the coordinate perpendicular to 



the surface. In Fs 1 , g and 7 are parameters which depend on temperature and the 
exchange couplings in the bulk (J) and at the surface (J s ) Q : 

(q-2)J s + J- k B T 



7 



k B (T c -T) 
J 
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T 



-1/2 



(4) 
(5) 



2i h k B {T c -T)' 

Here, g denotes the coordination number of a site, and is the bulk correlation 
length. [Our normalization of F[ip] in Eq. (j3J) implies that all lengths are measured 
in units of Further, the dimensionless surface field in Fg 1 is fti = — V(0). The 
one-sided derivative appears in Fs 1 due to the absence of neighboring atoms for 
z < 0. Similarly, the third term F$ 2 is the contribution from the surface S2 located 
at z — D, with h 2 = —V(D). For simplicity, we assume that J Sl = J s . 2 = J s , so that 
the parameter g is the same for both S\ and S2. 

The corresponding dynamical model is obtained as follows. In the bulk, the 
order parameter evolves according to the Cahn-Hilliard-Cook (CHC) equation for a 



conserved order parameter 



4l|: 



9 ,^ ^ 



-V • J(f, t) 

v • [Vii(f,t) + e(f,t)} 



(6) 



Here, J(r,t) is the current, and n(r,t) is the local chemical potential difference 
between A and B. Further, 9(f,t) is a random noise term, to be specified below. 
Using the free energy functional from Eq. © in Eq. (jHJ), we obtain 



J^(r,t) = V-<{V 



-if) + ip z 



-V 2 ^ + V(z) 



+ 6{f,t) 



< z < D. 



(7) 



We assume that the noise 9 is a Gaussian white noise, obeying the relations 

(9(r,t)) = , 
(e i (f',lf)e j (f",f)) = 2e5 ij 5(r'-r")5(t'-t") . 



(8) 
(9) 
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The dimensionless noise amplitude is a function of the temperature 

-2 







3 \T -• - ■ < 10) 
Eqs. ©-0 model the fact that the order parameter (the total concentration) is 
conserved (41 1. However, it is important to note that the surface value of the order 
parameter is not a conserved quantity. We assume a nonconserved relaxational 



kinetics (referred to as Model A 







4l|) for this quantity at Si. 
SF 



/n + 0V(P,O,*)+7 



dip 
dz 



where To sets the time scale of the nonconserved kinetics. Since the surface value of 
the order parameter relaxes much faster than the time scales of phase separation, 
we subsequently set r = 0. Finally, we observe that there is no current of material 
across S\. This is implemented via a no-flux boundary condition: 



Up,o,t) = -\j- z 



-ip + ip 3 --V 2 iP + V{z) 



o 



'121 



2 = 



The boundary conditions at z — D are implemented in a similar manner. For the 
sake of completeness, we present them here: 








h 2 + gip(p,D,t) -7 



dip 
dz 



z=D 



d_ 

dz 



-ip + ip 3 --V 2 ip + V(z) 



+ 



(13) 
(14) 



z=D 



Eqs. (J7 | )-()14 p constitute our model for phase separation in a film [39(. This model 
has been presented in the context of a film with flat parallel surfaces. However, the 
adaptation to an arbitrary geometry is obvious: boundary conditions like Eqs. (jlll)- 
(fT4^1 are implemented on all available surfaces. For example, Aichmayer et al. [42] 
have used the appropriate generalization of this model to study SDSD in a cylindrical 
geometry. 

We have undertaken a Langevin simulation of the above model, including the 
noise term, in order to study phase separation in both symmetric and antisymmet- 
ric films. We implemented an Euler-discretized version of Eqs. (|7])- (fT4^) on cubic 



lattices of size L 2 x D with L = 256 and D = 5, 10. The discretization mesh sizes 
were Ax = 1 and At = 0.02. We should stress that these mesh sizes are rather 
coarse and the resultant numerical solution does not closely shadow the "actual" 
solution of Eqs. ©-(JHJ). However, Oono and Puri j3] and Rogers et al. j^J have 
demonstrated that such discrete cell- dynamical system models capture the physics of 
the segregation process rather well. The boundary conditions in Eqs. (|lH ) - (fPH) were 
implemented at z = and z = D, respectively, while periodic boundary conditions 
were applied in the x- and y-directions. 

The nature of the surface potential V(z) and the parameters g,7,e determine 
the equilibrium phase diagram of the film. For the semi-infinite case, the phase 
diagram has been discussed in Refs. hJ 13, Q|. For the films considered here, we 
have determined the boundary between PW and CW phases both analytically and 
numerically. For the sake of brevity, we do not discuss these phase diagrams here, 
but point out that parameter values are chosen to study ordering to both PW and 
CW states in symmetric and antisymmetric films. The values of the parameters in 
the boundary conditions were chosen as g = —0.4, 7 = 0.4 ^J. The potentials 
V(z) which we considered will be specified at appropriate places in the subsequent 
discussion. The noise amplitude was fixed as e = 0.327, which corresponds to a 
quench with T ~ 0.38 T c from Eq. (fTT)|) . The presence of thermal fluctuations 
prevents the system from becoming stuck in metastable configurations. However, 
we should stress that the ordering dynamics is expected to be independent of noise 



in the asymptotic regime 



47]. 



The initial conditions for our simulations consisted of a homogeneous mixture of 
50% A and 50% B, i.e., ip(f, t = 0) = + small-amplitude fluctuations. This mimics 
the disordered high-temperature state for a mixture with critical composition. We 
will characterize the far-from-equilibrium dynamics of the quenched system via evo- 
lution snapshots, laterally averaged order parameter profiles, layer-wise correlation 
functions and length scales. All statistical quantities were obtained as averages over 
five independent runs. 
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III. KINETICS OF PHASE SEPARATION IN A SYMMETRIC FILM 



Let us first consider the case of a symmetric film with a power-law potential: 

V(z) = -V [(z + l)- 3 + (D + l-z)- 3 ] . (15) 

Recall that the exponent n = 3 corresponds to non-retarded van der Waals' interac- 
tions between a surface and a particle in d = 3. As the other parameters are fixed, 
an appropriate choice of Vq and D results in either PW or CW states in equilib- 
rium. We have ascertained the PW-CW phase boundary by studying the evolution 
of an initial condition which consists of A-rich and B-rich domains separated by an 
A-B interface along the ^-direction. The onset of the CW phase is signaled by the 
intrusion of a thin wetting layer [see Fig. E^b)] between the B-rich domain and the 
surfaces. We will consider the PW and CW cases separately. 



A. Partially Wet Surfaces 

In Fig. |2l we show evolution snapshots [part (a)] and (a;2;)-cross-sections [part 
(b)] for films with D = 5 (frames on left) and D — 10 (frames on right). The 
potentials were chosen with Vq = 0.325 for D = 5 and Vq — 0.11 for D = 10, 
which correspond to the PW state in equilibrium. (The critical value of Vo/a for 
the PW— >CW crossover diminishes with increase in D, and Vo/a — > 1 as-D^oo.) 
Note that a metastable layered structure forms at early times, since the kinetics of 
surface enrichment 4^ is much faster than the time scale of phase separation. On 
longer time scales, spinodal fluctuations break this layered structure and the system 
forms domains which coarsen in directions parallel to the surface. We stress that 
the layered state can be very long-lived, and may be misinterpreted as evidence for 
the formation of wetting layers in experiments. 

Many experimental probes (such as depth-profiling techniques) do not have any 
latera 



vs. z 



resolution, and yield only laterally averaged order parameter profiles ip av (z, t) 
I n our simulations, laterally averaged profiles are obtained by averaging 
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i/j(x, y, z, t) along the x- and ^-directions, and then further averaging over five inde- 
pendent runs. The depth profiles corresponding to the evolution in Fig.|2]are shown 
in Fig. El For bulk spinodal decomposition, the wave- vectors are randomly oriented 
and the averaging procedure yields ip av (z,t) ~ 0. For SDSD, however, the laterally 
averaged profiles are systematic at the surfaces since the boundary conditions result 
in spinodal waves with wave-vectors perpendicular to the surfaces. Let us focus on 
the case with D = 10 in Fig. E^b). The profile at time t = 10 shows the formation 
of two symmetric SDSD waves, which propagate towards the center of the film. The 
t = 100 profile shows the metastable layered state that has originated from these 
waves. This structure is also present at t = 1000, and may be misinterpreted as 
a CW equilibrium state which occurs for temperatures between the critical tem- 
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perature of the thin film and the critical temperature of the bulk system 
Finally, the spinodal fluctuations break this structure and the averaged profile at 
t = 20000 is almost flat. Since a weak surface field amplitude Vq = 0.11 was chosen 
in this case, there is only a slightly A- rich region [ip av (z,t) > 0] near the walls and, 
correspondingly, only a slightly A-poor region [%j) m {z, t) < 0] at the center. 

It is also interesting to study (:q/)-cross-sections of the evolution snapshots in 
Fig. 121 In Fig. we show the relevant cross-sections at z = 2 for D = 5 and z = 5 
for D = 10. For early times (t = 100) the central region is strongly depleted in A due 
to the formation of the layered structure. The resultant morphology corresponds 
to an off-critical composition with droplets of A in a matrix of B. At later times, 
t = 20000, the central region has almost equal amounts of A and B again. However, 
there is still a small depletion in A (see Fig. EJ), and hence the growth morphology 
still contains droplets of A. 



Let us next focus on the layer-wise correlation function, which is defined as 
CJp,z,t) =L~ 2 / da[(i:(a,z,t)ij(a + p,z,t))-(ij(a,z,t))(i}(a + p,z,t))} , (16) 



where the angular brackets denote an averaging over independent runs. Since the 
system is isotropic in the (x?/)-plane, C.. does not depend on the direction of p. We 
can define the z-dependent lateral length scale L (z, t) = L(z, t) from the half-decay 
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of C.. (p, z, t) 

C^p = L lZl t)= l -C^,z,t). (17) 

For convenience, we denote CJp,z,t) as C(p,t) in the following discussion. In 
Figs. ED^a),(b) we plot the scaled correlation functions C(p,t)/C(0,t) vs. p/L for 
D = 5 and D = 10, respectively. In bulk systems, the correlation function exhibits 
dynamical scaling, viz., C(f,t) = g(r/L), where g(x) is independent of time. This 
property indicates that the evolution morphology is statistically self-similar in time, 
and only the scale of the morphology changes. In this case, there is no dynamical 
scaling as the correlation functions correspond to qualitatively different morphologies 
(see Fig.|2J). Thus, for D = 10 and t = 1000, the layered structure has not yet broken 
up, while at t — 20000 lateral phase separation has occurred. Of course, dynamical 
scaling is recovered subsequent to the formation of well-formed laterally segregated 
domains (t > 10000 for both D = 5 and 10). 

Finally, we examine the time-dependence of the lateral domain size L(z,t) in 
Fig. El While the asymptotic growth is consistent with the expected Lifshitz-Slyozov 

(LS) growth law L(t) ~ t 1/3 |l] films 

, which describes bulk domain growth, the 
early-time dynamics is complicated. For the D = 10 case in Fig.Efb), the early-time 
data corresponds to the growth of bulk-like domains before the layered structure has 
formed. The spinodal fluctuations originate in the central region (z = 5) where the 
surface field is not felt, and propagate to the surface (z = 0). The break-up of the 
layered structure is characterized by the non-monotonic behavior of L(z,t) vs. t. 
As the fluctuations originate near the film center, the data set for z = is the last 
to become consistent with LS behavior. 



B. Completely Wet Surfaces 

Next, let us consider the case where the surfaces have a CW morphology in 
equilibrium. In Fig. we show evolution snapshots and (x2;)-cross-sections for the 
CW case. The corresponding potential strengths were Vq = 0.45 for D = 5, and Vq = 
0.275 for D = 10. Again, the system forms a metastable layered structure at early 
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times, which is broken up by spinodal fluctuations at later times. (For very strong 
surface fields, the layered structure actually corresponds to an equilibrium state and 
the corresponding pattern dynamics is uninteresting.) However, the difference in 
this case is that the B-rich regions are encapsulated by A [see Fig. ^b)], unlike the 
situation shown in Fig. |21 The asymptotic dynamics is then characterized by the 
coarsening of these encapsulated pancakes. The laterally averaged profiles (Fig. |HJ) 
again show that the initial layered structure is rather pronounced (compare with 
Fig. EJ). The depth profiles become softer at later times, but due to the pancakes 
being encapsulated by A, there remains a strong surface enrichment in A. If one 
looks at cross-sections taken parallel to the surfaces, analogous to Fig. 01 one finds 
a qualitatively similar behavior. Of course, the volume fraction of A in the central 
region is now smaller due to the higher degree of surface enrichment. 

FigureOlis a scaling plot of C(p, t)/C(0, t) vs. p/L, and is analogous to Fig. El We 
do not show data for the z = case here as the surface is always A-rich and does not 
exhibit interesting pattern dynamics. The behavior in the film center is qualitatively 
similar to the PW case, i.e., there is no dynamical scaling for the time range shown. 
This can be understood in the context of the evolution dynamics shown in Fig.[7fb). 
For the D = 10 case, the morphologies exhibit a crossover behavior from the layered 
state to the (asymptotic) pancake state for t = 100, 1000, 20000. For t > 20000, we 
expect to recover dynamical scaling. For the D = 5 case, the system is almost in 
its asymptotic state by t = 1000. Hence, the correlation functions for t = 1000 and 
t = 20000 show approximate scaling. 

In Fig.Hni we study the time-dependence of the parallel length scale L(z, t). The 
non-monotonic behavior again reflects the formation and break-up of a long-lived 
metastable layered structure. Notice that this state is far from equilibrium despite 
the fact that we do expect a CW morphology for these parameter values. The 
break-up of the layered structure gives rise to the growth of laterally segregated 
domains. Although we have followed the growth of L(z,t) over several decades in 
time, the expected asymptotic regime of LS growth is not observed over simulation 
time-scales. 
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IV. KINETICS OF PHASE SEPARATION IN AN ANTISYMMETRIC 
FILM 

We next consider the case of an antisymmetric film. The corresponding power-law 
potential is 

V(z) = -V [(z + iy 3 - (£> + 1 - z)~ 3 } , (18) 

so that V(D — z) = —V(z). In Figs, ^c)-(d), we have schematically shown the PW 
and CW states which arise for an antisymmetric film. In this case also, we have 
obtained the PW-CW phase boundary as a function of Vq and D. As before, we will 
consider and compare both PW and CW cases. 

A. Partially Wet Surfaces 

In Fig. HH we show evolution snapshots [part (a)] and (x2;)-cross-sections [part 
(b)] for films with D = 5 (frames on left) and D — 10 (frames on right). The 
potential strengths were Vq = 0.055 for D = 5 and Vq = 0.041 for D = 10, which 
correspond to a PW case in equilibrium. In the D = 5 case, we observe the formation 
of a layered state which breaks up into a coarsening domain structure. A similar 
evolution occurs in the D = 10 case, though the layered state (at t = 1000) is not so 
clean for these weak surface fields. However, it shows up more clearly in the laterally 
averaged profiles we present next. The domain cross-sections are trapezoidal with 
different contact angles at the lower surface (which prefers A) and the upper surface 
(which prefers B). 

The laterally averaged profiles (Fig. IT2|) confirm that a layered state, with a single 
interface, appears as a transient before the lateral domain growth sets in. For the 
D = 10 case, we see the interaction of two opposite SDSD waves (at t = 10, 100), 
resulting in the formation of the layered state (at t = 1000). In Fig. we examine 
the morphology in planes parallel to the surfaces. For the D = 5 case, we focus on a 
cross-section at z = 2. The t = 100 morphology corresponds to the layered state [see 
Fig. ll2f a)] and consists of droplets of B in a matrix of A. The t = 20000 morphology 
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corresponds to the laterally-segregated state and consists of bicontinuous domains. 
For the D = 10 case, we consider a cross-section at z — 5, which is precisely the 
film center. Since the film always has a near-critical composition at the center 
[see Fig. 112( b)], the segregation morphology is bicontinuous for both t = 100 and 
t = 20000 in this case. 

We have also examined the layer- wise correlation functions CJp,z,t) = C(p,t) 
f Fig. 114)). and the lengths L(z,t) which one can extract from them (Fig. ^J). Again, 
one finds marked deviations from scaling for C (p, t) , as expected due to the transient 
formation of layered structures. Of course, scaling is recovered in the asymptotic 
regime, which is characterized by the coarsening of trapezoidal domains. In the 
D = 5 case, the length-scale data shows that the asymptotic behavior is consistent 
with the LS law, L{t) ~ t 1 / 3 . A similar behavior is seen for the film with D = 10, 
though the non-monotonic behavior is less pronounced in this case. This emphasizes 
that one has to be careful with the interpretation of growth phenomena in confined 
geometries, and rather complete information on the structural evolution of a system 
is mandatory for establishing a clear picture. 

B. Completely Wet Surfaces 

Our last set of numerical results corresponds to the CW state for an antisym- 
metric film: in this layered state with a single interface is the equilibrium 
state [see Fig. GJd)], and no lateral segregation should occur! In Fig. HHJ we show 
evolution snapshots and (a;z)-cross-sections for D = 5 with Vo = 0.25 (frames on 
left) and D = 10 with Vo = 0.2 (frames on right). There is seen to be some lateral 
inhomogeneity in the early stages. Starting from a random initial state, small A- 
rich and B-rich domains are formed first. It takes time for the interfaces between 
these small domains to annihilate by diffusion and coalescence, until a single domain 
wall parallel to the surface is left. The corresponding laterally averaged profiles are 
shown in Fig. El In particular, we draw the reader's attention to Fig. ITTT b). which 
shows the formation and collision of two opposite SDSD waves originating from Si 
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and S2 - see profiles for t = 10 and t = 100. Figure fTSl shows the evolution pictures 
at the film center (z = 5) for the D = 10 case. (The corresponding pictures at z = 2 
or z = 3 for the D = 5 case merely show a uniform state.) The length scale at 
z = 5 for £) = 10 (Fig. IT^j) grows uniformly for t > 1000, and the time-dependence 
is consistent with the LS law. 



V. SUMMARY AND DISCUSSION 

Let us conclude this paper with a summary and discussion of the results. We 
have studied the diffusion-driven phase separation of an AB mixture confined in a 
film. The film has two parallel surfaces Si and S2, which are separated by a dis- 
tance D ~ O(10£b), where is the bulk correlation length. We have considered (a) 
symmetric films, where Si and S2 have an identical attraction for the A-component; 
and (b) antisymmetric films, where S\ and S2 attract (with equal strength) the 
A-component and B-component, respectively. Both cases are of considerable exper- 
imental relevance. 

The equilibrium segregated state can be either partially wet (PW) or completely 
wet (CW), depending on the nature of the surface potentials. Further types of 
mixed-phase states in thin films occur only for restricted ranges of parameters 3J, |35 , 
[36^ , and are not considered in the present paper. We have clarified the typical growth 
scenario in both symmetric and antisymmetric films. In both PW and CW cases, the 
surfaces give rise to surface-directed spinodal decomposition (SDSD) waves, which 
propagate towards the film center. The interaction of these SDSD waves leads to the 
formation of a layered state. This state is metastable for the PW case, and is broken 
up by spinodal fluctuations. (However, the metastable state may have a very long 
lifetime, and could be of considerable experimental significance.) The subsequent 
evolution of the mixture is characterized by the lateral coarsening of pancake-like 
domains. (For antisymmetric films, these domains are trapezoidal because of the 
different contact angles at Si and 52.) In the later stages, we expect that this 
coarsening is governed by the Lifshitz-Slyozov (LS) growth law L{t) ~ t 1 / 3 , but one 



15 



often encounters slow transients and even a non-monotonic variation of L(t) with 
time, before the LS regime sets in. For the CW case in symmetric films, the initial 
layered state is again metastable for moderate surface fields, and breaks up into 
encapsulated pancakes which coarsen in the lateral direction. For the CW case in 
antisymmetric films, the equilibrium state is a layered state with a single interface 
parallel to the surfaces. The system can relax to this state rather rapidly. 

At this stage, one may ask what happens if hydrodynamic effects are incorporated 
into the above discussion. It is well-known that bulk fluid mixtures 

exhibit more complicated segregation kinetics than solid mixtures. In the initial 
stages, growth is diffusive and is governed by the LS growth law. However, at later 
times, hydrodynamic effects become relevant and domain growth is facilitated by 
advective transport along interfaces. The corresponding growth laws are L(t) ~ t 
in the viscous hydrodynamic regime and L(t) ~ t 2//3 in the inertial hydrody- 
namic regime [50| • The effects of surfaces on phase-separating binary fluids can be 
studied at various levels of description. At the coarse-grained level, an appropriate 
model is Model H at a surface 0|. This consists of the coupled dynamics of 
an order-parameter field and a velocity field, with appropriate boundary conditions 
at the surfaces. Alternatively, one can study mesoscale models consisting, e.g., of 
evolution equations for the configuration probability distribution j^j]. Finally, at 
the microscopic level, one can undertake molecular dynamics (MD) simulations of 



binary fluids 



53| in a confined geometry. We have undertaken such MD simulations, 



and will report the results in a forthcoming publication. 

More generally, we emphasize that there are many intriguing aspects of phase 
separation in confined geometries, and a number of possible directions for further 
investigation. For example, the present study focused on mixtures with critical com- 
position. In a semi-infinite geometry, Puri and Binder j3| have demonstrated that 
novel features arise when off-critical compositions are considered, e.g., the wetting 
layer grows faster when the wetting component is a minority phase rather than 
a majority phase. It would also be interesting to study the phase separation of 
off-critical mixtures in the present context of confined geometries. Another inter- 
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esting complication arises if there is frozen-in disorder at the surfaces, which affects 
the wetting behavior. Finally, it is also relevant to study SDSD in more complex 
confined geometries than those studied here, e.g., wedges, patterned surfaces, etc. 
Though the modeling of these problems is straightforward, we expect that they will 
give rise to novel physical phenomena. There remain many issues to be addressed 
in this area, and we urge experimentalists to undertake fresh experiments in these 
directions. 
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Figures and Figure Captions 



a) PW 



1 


3 


1 


L M 


b) CW 


* 2R 


1 


A 


I 


B 


c)PW 


« 2R 




1 




\ m 


d) CW 


2R ► 


B 




Aj 





D 



D 



D 



D 



FIG. 1: Schematic description of phase-separated states in thin films of volume L x L x D. 
We show cross-sections of the film in the (xz)-plane. Cases (a)-(b) refer to thin films where 
both surfaces prefer A, and the surface potential is symmetric, V(D — z) = V(z). Cases 
(c)-(d) refer to thin films where the lower surface prefers the A-rich phase and the upper 
surface prefers the B-rich phase. The corresponding surface potential is antisymmetric, 
V(D — z) = —V(z). The A-rich domains are marked black, and the B-rich domains 
are unmarked. For both symmetric and antisymmetric films, partially wet (PW) and 
completely wet (CW) morphologies emerge in the limit D — > oo. Further types of phase- 
separated states exist in thin films at off-critical composition [34L 13.1 13n| . but we will not 
consider these here. Note that the thickness d of the encapsulated B-rich domain in (b) 
differs from D only by corrections which increase slower than linearly with D. 
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FIG. 2: (a) Evolution snapshots at times t = 100 and t = 20000 for a critical binary 
mixture in a symmetric film with PW morphology. The system size was L 2 x D with 
L = 256 and D = 5 (frames on left) and D = 10 (frames on right). The A-rich regions 
are colored black, and the B-rich regions are colored white (light blue online), (b) Per- 
pendicular cross-sections of the snapshots in (a) at y = L/2 in the (xz)-plane. The A-rich 
regions are marked black, and the B-rich regions are unmarked. 
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(a) D=5 (PW) 




z 



FIG. 3: Laterally averaged profiles for the evolution depicted in Fig.|2]at the dimensionless 
times t = 10, 100, 1000, 20000, for (a) D = 5, and (b) D = 10. The symbols denote the 
same times in both figures. 
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FIG. 4: Cross-sections of the evolution snapshots in Fig. |2fa). The cross-sections are 
taken parallel to the surfaces at (a) z = 2 for D = 5, and (b) z = 5 for I? = 10. 
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FIG. 5: Scaling plot of layer- wise correlation functions, denned in Eq. (|lfjjl. for the evolu- 
tion depicted in Fig. We plot data for C(p,t)/C(0,t) vs. p/L for three different times, 
t = 100, 1000, 20000. We present data for (a) D = 5 at z = (wall) and z = 2 (center); 
and (b) D = 10 at z = (wall) and z = 5 (center). The layer-wise length-scale L(z,t) 
is defined as the distance over which C(p,t) has decayed to 1/2 its maximum value (at 
p = 0). 
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FIG. 6: Time-dependence of the layer-wise length scale for the evolution depicted in 
Fig. [2j We plot L{z, t) vs. i on a log-log scale for various values of z and (a) D = 5, 
and (b) D = 10. The lines of slope 1/3 denote the Lifshitz-Slyozov (LS) growth law, 
L(t) oc i 1 / 3 , which characterizes diffusion-driven phase separation in the bulk. 
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FIG. 7: Analogous to Fig. but for a symmetric film with a CW morphology. The 
parameter values are just above the PW-CW boundary for the D = 10 case. Therefore, 
the surface does not have a clean coating of the A-rich phase. 
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FIG. 8: Laterally averaged profiles for the evolution depicted in Fig. dat the dimensionless 
times t = 10, 100, 1000, 20000, for (a) D = 5, and (b) D = 10. 
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FIG. 9: Scaling plot of layer- wise correlation functions for the evolution depicted in Fig.0 
We plot C{p,t)/C(0,t) vs. p/L for t = 100, 1000,20000. We present data for (a) D = 5 
at z = 2 (center); and (b) D = 10 at z = 5 (center). 



29 



(a) D=5 (CW) 



10: 



N 




10: 



N 



(b) D=10(CW) 

i 1 1 1 1 1 mi i 1 1 1 1 1 mi i 1 1 1 1 1 mi i 1 1 1 1 1 mi i 



o z=3 
□ z=4 
o z=5 



1/3 



o 



°oooo ooocP 0000 JJ 



,9 



i i i i m mi i i i 1 1 1 hi i i i 1 1 1 hi i i i 1 1 1 hi r 

f) 1 2 3 4 

10 10 10 10 10 

t 



FIG. 10: Time-dependence of the layer-wise length scale for the evolution depicted in 
Fig. [7J We plot L(z, t) vs. t on a log-log scale for various values of z and (a) D = 5, and 
(b) D = 10. 
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FIG. 11: Analogous to Fig. but for an antisymmetric film with a PW morphology. 
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FIG. 12: Laterally averaged profiles for the evolution depicted in Fig. ^2 at the dimen- 
sionless times t = 10, 100, 1000, 20000, for (a) D = 5, and (b) D = 10. 



32 



(a) D=5 (PW) 

t=100 



(b) D=1 (PW) 

t=100 





*- *' k< 




1 




x V . ■ 


























f ■ 




ft 






- 










* 




-., 


* 1 * 


* t 


* > 


1 

■■ -V 




















if 




' . 1 y 




1 * 












• 






-* J w 


* * 




+ 


.+ 




^ ■ 


■ V 

If 




.* '4 * 


\ 

■ * A 








33 



(a)D=5(PW) (b)D=10(PW) 
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FIG. 14: Scaling plot of layer-wise correlation functions for the evolution depicted in 
Fig.HH We plot C(p,t)/C(0,t) vs. p/L for t = 100, 1000,20000. We present data for (a) 
D = 5 at z = (wall) and z = 2 (center); and (b) D = 10 at z = (wall) and z = 5 
(center). 
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FIG. 15: Time-dependence of the layer-wise length scale for the evolution depicted in 
Fig.HU We plot L(z,t) vs. t on a log-log scale for various values of z and (a) D = 5, and 
(b) D = 10. 
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FIG. 16: Analogous to Fig. but for an antisymmetric film with a CW morphology. 
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FIG. 17: Laterally averaged profiles for the evolution depicted in Fig. Ilfil for (a) D = 5, 
and (b) D = 10. 
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FIG. 18: Cross-sections of the evolution snapshots for D = 10 in Fig. 116( a) . The cross- 
section is taken parallel to the surfaces at z = 5. 
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FIG. 19: Time-dependence of the layer-wise length scale for the D = 10 evolution depicted 
in Fig. El We plot L(z, t) vs. t on a log-log scale for z = 5. 



39 



